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Abstract 

We calculate QNMs of the scalar hairy black hole in the AdS back- 
ground using Horowitz-Hubeny method for the potential that is not 
known in analytical form. For some black hole parameters we found 
pure imaginary frequencies. Increasing of the scalar field mass does 
not cause the imaginary part to vanish, it reaches some minimum and 
then increases, thus in the case under consideration the infinitely long 
living modes (quasi-resonances) do not appear. 

1 Introduction 

There are three stages of the perturbation evolution near a black hole. The 
initial black hole outburst from the perturbation source takes place at early 
times of the evolution. Then at late times the damping (quasi-normal) oscil- 
lations appear, finally suppressed by power-law or exponential tails at asymp- 
totically late times. The quasi-normal modes (QNMs) can be characterized 
by a set of complex frequencies: real part is the actual oscillation frequency 
while imaginary part describe the damping rate of the oscillation. The QNMs 
do not depend on initial perturbations being, thereby, an important charac- 
teristic of the black hole. 

For this reason, QNMs of black holes in asymptotically flat space-time 
have been extensively studied for more then thirty years (see reviews [1J). 
Yet the quasi-normal spectrum of a black hole is crucially dependent on a 
cosmological background which a black hole is immersed in. QNMs of black 
holes in de Sitter background were considered recent years in jlj |3] . 

Despite we live in the de Sitter universe, black holes in the 
anti-de Sitter(AdS) background attracted considerable interest due to the 
famous AdS/CFT correspondence |3 4 : it is shown that large black holes 



in AdS space-time correspond to approximately thermal states in the dual 
conformal field theory(CFT). According to this correspondence, QNMs of 
a black hole should coincide with poles of retarded Green function of the 
perturbation in M = 4 SU(N) super- Yang-Mills theory for large N and were 
extensively studied recent years jU E] . 

First of all, let us briefly review what is known about perturbations of a 
scalar field interacting with a black hole. Scalar field is the simplest model, 
frequently used, when we neglect spin influence [HJ E]- Also scalar field is 
used to explain dark matter problem (see e. g. [ID])- In fact almost all 
papers which considered scalar field perturbations and corresponding quasi- 
normal modes were limited by minimal coupling of the scalar field with a 
black hole, i. e. by consideration of test scalar field propagating in the black 
hole background fTJ ^] . Further extensions within minimal coupling were 
made to consider massive scalar field El ■ 

Very few papers beyond the minimal coupling are works where perturba- 
tions of a black hole with specific form of scalar field called dilaton were con- 
sidered |15( ITS] , and also perturbations of Bocharova-Bronnikov-M elnikov- 
Bekenstein (BBMB) black hole [T2j studied in |T8|. Note however, that 
BBMB black hole is topologically equivalent to extremal Reissner-Nordstrom 
solution and therefore does not provide any unexpected features [THj . 

Since the paper of J. Bekenstein |19j . it is well-known that black hole can 
not have scalar hair within minimal coupling. Yet as recently have been found 
in (221 121] , there is a possibility of dressing a four-dimensional black hole in 
AdS space-time with a non-minimally coupled classical scalar field. In 25 
the results were extended for higher dimensional configurations. Moreover, 
fortunately the found solution is stable, thereby allowing us to study its 
quasi-normal modes. Thus one motivation to undergo this study is to find 
out what will happen with the ringing behavior of a black hole with scalar 
hairs, Another aim is rather technical: for the above black hole perturbations 
the problem is reduced to wave-like equations but with effective potential 
which is not known in analytical form, but expressed in terms of unknown 
functions obeying some differential equations. 

This is the case of many physically relevant situations such as dirty black 
holes (i. e. with some distribution of matter around them), black holes with 
massive dilaton, evaporating black holes and other objects determined by a 
set of equations that can not be solved analytically. We hope therefore that 
the way we are considering quasi-normal spectrum here may be of some use 
in those cases. 
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The other interesting issue is related to arbitrary long living modes (in lin- 
ear approximation) which were called quasi-resonances. They appear when 
perturbing massive scalar |2DJ or vector t 2T] fields near a black hole of a 
specific mass. In [22] it was shown, that quasi-resonances are permitted for 
black holes immersed in asymptotically flat background and are not permit- 
ted for asymptotically de Sitter space-times. We show that anti-de Sitter 
boundary conditions do not permit quasi-resonances too (see the appendix). 
Our numerical investigations confirm this statement. 

This paper is organized as follows. In section |2] we review briefly [23121] 
to introduce the hairy black hole and the perturbation equation near it. In 
section |H] we describe the method we used for our calculations. And finally, 
in section H] we make some comments about the results. 

2 Scalar hairy black hole 

We consider spherically symmetric solution 

ds 2 = N(r)e 2S{r) dt 2 - iV(r)-Mr 2 - r 2 (d6 2 + sin 6 2 dip 2 ) (1) 

of the action, which describes a self-interacting scalar field <fi with non- 
minimal coupling to gravity: 



S = d 4 x y/—g 



\ (R - 2A) - ~ (V0) 2 - \ m 2 ~ V{4>) 



(2) 



where R is the Ricci scalar, A is the cosmological constant, £ is the coupling 
constant, V{<p) = — - — is the massive scalar field self-interaction potential 

and (V0) 2 = V M 0V^0. 

In [21] stability of the non-trivial solution has been proven for < £ < 
3/16, 1 — £0 2 > 0, fi = 0. In 23 j some arguments for its stability for 
< /i 2 < — 4£A are given and the stability is proven for the particular 
case of £ = 1/6. In the present paper we consider the above range of the 
parameters and our results support stability of the solution within linear 
perturbatuon approach. 

It is convenient for our purposes to use conformal transformation in this 
region to simplify the equations of motion: 

g^v = (l - ^4> 2 )g^- (3) 
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After the transformation the action takes the form: 



S = |^- v ^Q(i?-2A)-i(v$) 2 -f/(«l>)), (4) 



where we define 
$ = 



1 - m + 6£ 2 2 _ V(<f>) + A^ 2 (2 - 

\ (i - m 2 ' { } ~ o^f ' { } 

The metric takes the following form 

ds 2 = N(f)e 2 ~ m dt 2 - N{f)- 1 df 2 - f 2 (dO 2 + sin # 2 dip 2 ) , (6) 

where 



r = ( 1 — f d) 2 ) 2 r, 



2/ - ^2r 2 



AT = AT (l - £0 2 - £ r 00') (l 
Ne 25 = Ne 25 (l - £0 2 ) . 

Varying the action one can find the equations 

d(fN) k o o ( N (d§\ 2 

)\ (8) 

(n— — —\ — — 

\ df dr f J df d& 

In order to obtain the explicit form of © we have to solve the equations 
with some boundary conditions that are associated with the black hole pa- 
rameters. We are using the following parameters. The event horizon radius 
fh can be chosen arbitrary in order to fix length scale and defines boundary 
condition for N: N(fh) — 0. We measure all dimensional units in units of f^, 
thus we choose f^ — l. S(fh) can be also chosen arbitrary to fix time scale. 
In order to introduce the boundary conditions in the same point we choose 
^(^h) = as distinct from Winstanley's S(oo) = 0. And the last parameter 
4>{fh) — Q can be associated with the scalar charge of the black hole. 



d6 




df 


2 V ~dJ 











df 2 
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Spherically symmetric perturbation equation for exp(— iujt)%p(r) = r<5$ 
takes the Schrodinger wave-like form: 



d*if> 



+ (u 2 - tf)if> = 0, (10) 



dr* 2 

where the perturbation potential U is given by 

Ne 2S {, , TT > o \ ( o{d<!>\ 2 \ - ~dUd$ 2 d 2 U 



r* is tortoise coordinate given by 

dr 



and similar to ordinal anti-de Sitter case: 



r = rh -<=>- r = oo, 
f = oo <^=^> f* = 0. 



:u) 



e'iV, (12) 



This perturbation corresponds to infinitesimal changing of the black hole 
mass and is similar to the zero-multipole scalar perturbations near a hairless 
black hole. The crucial difference is that the perturbation of a scalar field 
due to the hair gives a first-order correction to the metric. 



3 QNMs calculation 

Quasi-normal modes can be defined as eigenvalues of (jl(Jj) with the boundary 
conditions which require ingoing wave at the event horizon and zero at spatial 
infinity. These boundary conditions are natural and mean that the black 
hole does not emit radiation and our solution is not divergent. These specific 
boundary conditions allow us to fix the one of the two independent solutions 
in both points and calculate eigenvalues by comparing series near each of 
them. This technique is the base for two main precise methods of numerical 
calculation of QNMs: the continued fraction method 8\ and the Horowitz- 
Hubeny method [2J. They both were extensively used to find eigenfrequencies 
of perturbation equations with potential given analytically as polynomial 
fraction. Fortunately Horowitz- Hubeny method is possible to be adopted for 
equations with unknown potential. 
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Indeed if the potential U would be polynomial fraction we could eliminate 
singularity at the horizon introducing new function 



y(z) = z^(z), Z = 1 -J> ( 13 ) 

where k is chosen to make y(z) regular at the horizon z — 0. The require- 
ment for y(z) to be regular at the horizon fixes that independent solution 
which satisfies the horizon boundary condition. To use the other boundary 
condition one needs to find an equation for y(z) in the form: 

s(z)z 2 (j)"(z) + t(z)z(f)'(z) + u{z)(j){z) = 0, (14) 

where s(z), t(z) and u(z) are polynomial functions. Using explicit form of 
the polynoms: 

M s Ah M u 

S(z) = S nZ n , t( Z ) = J2 tnZ n , u(z) = ^ «n Z ", 

n=0 n=0 n=l 

it is easy to find series expansion for y(z) = 1 + ' s ^ j y n z n (we suppose that 

n 

there is no singularities in the circle \z\ < 1, if singularities appear one needs 
to use another definition of z to avoid them): 

^ n(n-l)s + nt 

and find uj by minimizing |y(l)| restricted by some large number N of ex- 
pansion terms in ()15|) . 

It is clear that if s(z), t(z) and u(z) are series and s(l), t(l) and u(l) 
converge quickly enough y n is still possible to calculate because it depends 
insignificantly on higher terms of the series. 

Thus we use the following procedure of the QNMs calculation. We 
introduce g(f) = N(f)/f 2 that is finite in the region outside black hole 
fh < f < oo. Using (J1(J|) we express s(z), t(z) and u(z) in terms of g{z), 
cf)(z) and 5(z) and find that k = l/g'(0). 

From ((71 |S1 EH) we calculate for given black hole parameters the first M ^> 1 
coefficients in the series expansions for g(z), <ft(z) and 5(z). Then for each 
given uj we find the first M coefficients of s(z), t(z) and u(z) and using (fTo'j) 
the first N ^> M terms of y(z). And, finally, we minimize \y(l)\ with respect 
to to. 
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We increase M and N for every given M until the value of u stops chang- 
ing within required precision. If we find finite M and N for a given precision 
of lu, we believe that the procedure converges and the u is eigenvalue of (fTU|). 
To control convergence of the procedure we used the simplest minimization 
method: we start from some initial value of u and calculate y(l) for this 
value, then we change uj by small value and calculate y(l) for it doing so 
while its modulus is decreasing, after that we decrease the small value and 
repeat the procedure until the required precision is reached. 

Obviously, the series s(z), t(z) and u(z) are definite up to an arbitrary 
finite function. Since we require quick convergence of them the selection of 
such function can greatly decrease M and N and therefore the calculation 
time. We used (1 — z) k factor for this purpose. (The complete procedure 
of the QNMs searching written in MATHEMATICA is available from the 
author upon request.) 

4 Discussion 

As were noticed above the method is appropriate only if equation (jl4j) has no 
singular points in the circle \z\ < 1. As long as functions s(z), t(z) and u(z) 
are known only as series we can not be sure that it is true within our approach. 
We believe that it is true because our results are stable against small changing 
of any of free parameters. If the equation had singularities in the unit circle 
we would obtain random results that do not approximate a smooth curve. 
Unfortunately we have no results to compare with ours because the hairy 
black hole is found as non-trivial solution of non-minimally coupled scalar 
field to gravity and has no studied limit to compare. Therefore calculation 
with an alternative method would be welcome to check the obtained results. 

Taking into account these argues we restricted our calculations to fun- 
damental modes 1 though the method allows to find higher ones that are 
presented for particular values of the parameters (see figure UJ). One can see 
that higher overtones tend to equidistant spacing similar to hairless black 
hole spectrum behavior in AdS [T3]. Yet fundamental modes give us a qual- 
itative picture of quasi-normal spectrum depending on the free parameters: 
£, Q, A and ji (we measure all quantities in ?h units). 

1 Fundamental QNM is the eigenfrequency with the lowest absolute value of its imagi- 
nary part. 



7 



Imaginary part demonstrates the most predictable behavior. It changes 
monotonously with respect to A, Q and £. Scalar field mass cause imaginary 
part to increase but not to vanish. It reaches a maximum for some particular 
mass and then decrease very rapidly (see figure EJ). Thus no infinitely long 
living oscillations appear. This fact is not unexpected and seems to be general 
for all black holes in the asymptotically anti-de Sitter background. 

Real part has more complicated behavior and depends significantly on 
the parameters Q and A (see figures El and EJ) • Other parameters £ and \i 
change the real part within the comparable small bounds (figures H] and • 
For relatively large values of £, A, Q purely imaginary frequencies appear. 

Also we did not find any growing modes (with positive imaginary part) 
that supports stability of such hairy black holes. 



Appendix: No quasi-resonance existence proof 

In j22] was shown that the existence of purely real quasi-normal modes (quasi- 
resonances) of massive scalar field near Schwarzschild black hole is due to 
specific boundary conditions at spatial infinity for the case of asymptotically 
flat space-time. It was explained also why there is no quasi-resonances of 
massive scalar field in the asymptotically de Sitter background if the potential 
is zero at both the horizons. Though the potential kind and the boundary 
conditions for asymptotically anti-de Sitter space-time are different one can 
show similarly that no quasi-resonances can appear in this background too. 

Indeed, suppose u is real and consider an imaginary part of the integral 
of (HDD: 



f*=0 



r*=—oo 



dip 
df* 



dr 



f*=0 



+ 



r"=—oo 



f*=0 



r'=— oo 



dip 



dr* 



+ {u 2 -U)\iP\ 2 \df*=0. (16) 



Because U is real 



= Im ip 



* dip 
df* 



=o r*=a 
+ 



r"=— oo 



dip 



dr* 



+ (u 2 -U)\iP\ 2 \ df*\ = (17) 
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Since we require Dirichlet boundary conditions at spatial infinity and by 

e*N 

1 ~ A eff, 



definition lim 



dip 
df* 



ip*e°N 



f*=0 



dr 



is finite. Indeed convergence of (|13|) at spatial infinity implies that 



V = z~^y{z) = 1 



-iu!K °0 



n=0 



n oo 



n=l 



;is) 



where C n can be expressed in terms of y n and thus e iV 
Therefore 



dr 



-CiA 



Im I ip*^r- 
dr* 



f*=0 



-Im f 



(19) 



that is impossible because ?/o = leads to trivial solution. Thus we prove 
that there can be no quasi-resonances. 
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Figure 1: First five quasi-normal frequencies for A = —0.1, Q = 1, 
£ = 0.1, /i = 0. They approach equidistant spacing as well as ordinary blach 
hole ones in AdS space-time. 



■ Im a) 



[23] E. Winstanley, Found. Phys. 33(2003) 111-43; ^^^^70 205092 

[24] E. Winsta nley Class. Quant. Grav. 22 (2005) 2233-2248; 
|gr-qc/0501096| 

[25] E. Radu, E. Winstanley Phys.Rev. D72 (2005) 024017; |gr5qc70503095] 



13 




(A = -0.05, f = 0.1, (i = 0) 




(A = -0.1, f = 0.1,(1 = 0) 

Figure 2: Dependance of the real and imaginary part of QN fre- 
quency on Q. The imaginary part quickly decreases with increasing of the 
"scalar charge" Q. The real part reaches its maximum and then quickly falls 
down to zero. It vanishes for some non-critical charge Q < C~ 1/2 that de- 
pends on other parameters: A, £, (i. For Q > Qq the real part remains zero 
and the frequency is purely imaginary (within numerical precision). 
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(Q = l, e = 0.1,/i = 0) 



Re w 




(g = 2.2, £ = 0.1, ii = 0) 

Figure 3: Dependance of the real and imaginary part of QN fre- 
quency on A. The real part increases for relatively small "scalar charges" 
Q, but for large Q = 2.2 it also reaches maximum and then falls down and 
vanishes A m —0.292. For larger values of A the real part remains zero and 
purely imaginary frequencies exist. The imaginary part decreases almost lin- 
early with A. You can see a small deflection from linearity on the bottom 
graphic. This deflection appears in the region where the real part decreases. 
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Im to 




Figure 4: Dependance of the real and imaginary part of QN fre- 
quency on £ < jq. Dependance of the quasi-normal frequency on £ is 
different for different values of A and Q. The imaginary part decreases with 
increasing of £. The quickness of such decreasing depends on A and Q. For 
some A and Q (as presented on the bottom graphic) the real part can fall 
down to zero and purely imaginary frequencies appear for larger £. 
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(A = -0.05, Q = 2.2, £ = 0.1) 

Figure 5: Dependance of the real and imaginary part of QN fre- 
quency on pi. We consider only < fi 2 < — 4£A. The real part changes 
within comparatively small range. The imaginary part reaches its maximum 
remaining negative. For this value of \x the oscillations have the longest 
lifetime, but they still damp and we do not observe quasi-resonanses. For 
larger scalar field masses the imaginary part of frequency decreases and the 
damping rate is higher. 
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